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We present a numerical study of an Ising spin glass with hierarchical interactions—the hierarchical 
Edwards-Anderson model with an external magnetic field (HEA). We study the model with Monte 
Carlo (MC) simulations in the mean-field (MF) and non-mean-field (NMF) regions corresponding 
to d > 4 and d < 4 for the d-dimensional ferromagnetic Ising model respectively. We compare 
the MC results with those of a renormalization-group (RG) study where the critical fixed point is 
treated as a perturbation of the MF one, along the same lines as in the e-expansion for the Ising 
model. The MC and the RG method agree in the MF region, predicting the existence of a transition 
and compatible values of the critical exponents. Conversely, the two approaches markedly disagree 
in the NMF case, where the MC data indicates a transition, while the RG analysis predicts that 
no perturbative critical fixed point exists. Also, the MC estimate of the critical exponent v in the 
NMF region is about twice as large as its classical value, even if the analog of the system dimension 
is within only ~ 2% from its upper-critical-dimension value. Taken together, these results indicate 
that the transition in the NMF region is governed by strong non-perturbative effects. 
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The renormalization group (RG) is a powerful method 
to tackle the complexity of a variety of thermodynami¬ 
cal systems close to their critical point: with the RG, an 
extensive number of microscopic variables is reduced to 
a few coarsened degrees of freedom whose behavior can 
be completely characterized in the thermodynamic limit. 
The RG theory for homogeneous systems—such as the d- 
dimensional Ising model on a hypercubic lattice—is based 
on a perturbative expansion around mean-field (MF) the¬ 
ory. Corrections to the MF fixed point (FP) are described 
in terms of a power series in e = 4 — d [l| , and the pre¬ 
dictions of this expansion in three dimensions have been 
shown to be in excellent agreement with both numerical 
simulations [2j) and experimental measurements from a 
wide variety of systems—binary fluids, superfluid helium, 
and ferromagnets Q. For inhomogeneous, disordered sys¬ 
tems characterized by a complex energy landscape such 
as Ising spin glasses [I|, the formulation of a RG theory 
is more involved: perturbative RG studies based on the 
e-expansion 0,0 are limited by the intricate diagram¬ 
matic structure of the underlying field theory, and the 
markedly divergent behavior of the e-expansion prevents 
any physical predictions below the upper critical dimen¬ 
sion [7]. 

In the last few decades, perturbative RG approaches 
have been extensively applied to a disordered system 
of particular relevance in condensed-matter physics: an 
Ising spin-glass in the presence of an external magnetic 
field [8]. Indeed, the existence of a spin-glass transition in 
a field has been raising particular interest [M3 because 
evidence in favor or against this transition would provide 
new insights into the structure of the low-temperature 
equilibrium states of the system [ 1 .'it . The spin-glass tran¬ 
sition in a field has been extensively studied with RG 
methods by examining a tentative set of critical FPs 0 - 
fl6j , although the existence of such a transition below the 
upper critical dimension is still under debate m, H3i- in 
this regard, a recent RG study [l8} has shown for the first 


time that the complete set of perturbative FPs in a field 
can be determined for a spin glass built on a hierarchi¬ 
cal lattice 0| : besides the well-known trivial FP in the 
MF region— corresponding to d > 4 for a ferromagnetic 
Ising model in d dimensions—this analysis demonstrated 
by means of an e-expansion that there is no perturbative 
FP in the non-mean-field (NMF) region 0|. The latter 
result provides a natural starting point for testing non- 
perturbative effects in spin glasses in a field: given that 
the RG analysis determined the complete set of pertur¬ 
bative FPs, if another FP in the NMF region were found 
with a non-perturbative approach, then this FP would 
necessarily lie outside the perturbative regime. 

In this paper, we investigate these non-perturbative ef¬ 
fects by studying with Monte Carlo (MC) simulations the 
hierarchical Edwards-Anderson model in a magnetic field 
(HEA)—a pairwise spin-glass model with binary cou¬ 
plings and magnetic fields where spin couplings are dis¬ 
posed in a hierarchical way. The HEA and the hierarchi¬ 
cal model in the RG study above [l8| share the same dis¬ 
tance dependence of the coupling strength between spins: 
this common feature between the two models allows us 
to compare the predictions of MC simulations with those 
from the RG approach. Specifically, we study the exis¬ 
tence of a phase transition by analyzing the temperature 
dependence of the correlation length and spin-glass sus¬ 
ceptibility, and we characterize the value of the critical 
exponent v [![]. In the MF region, the MC data and the 
RG approach both predict that there is a transition in a 
field, and that this transition is characterized by a classi¬ 
cal value of the critical exponent. Conversely, in the NMF 
region the two approaches exhibit a marked discrepancy: 
The RG method shows that there is no perturbative FP, 
while the MC data hints that there is a transition. In ad¬ 
dition, the MC value of the critical exponent v differs by 
nearly a factor of two from its classical RG value even if 
the analog of the system dimension is within only ~ 2% 
from its upper-critical-dimension value. Overall, these 
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findings suggest that the phase transition in the NMF re¬ 
gion is associated with a FP that cannot be considered as 
a perturbation of the MF FP, thus hinting at the presence 
of strong non-perturbative effects. 


characterized by vanishing order-parameter fluctuations 
for 1/2 < a < 2/3, and a NMF behavior in the region 
2/3 < cr < 1 where order-parameter fluctuations arise 


I 


The rest of the paper is structured as follows: In Sec¬ 
tion [T] we introduce the HEA, and in Section HAl we dis¬ 
cuss how a suitable correlation length based on the hi¬ 
erarchical distance can be extracted from the correlation 
functions. In Section HBl we present the MC data for the 
correlation length and overlap probability distribution, 
and in Section ED we discuss the numerical evidence for 
non-perturbative effects. Finally, in Section [II] we discuss 
the physical interpretation of our results and we lay out 
some topics of future studies. 


We will now study the existence of a phase transition 
related to the divergence of the spin-glass susceptibility 
in the HEA: to do so, we will introduce the correlation 
length associated with the long-wavelength modes of the 
spin-glass susceptibility and we will study it as a function 
of temperature for different system sizes. 


A. Correlation length 


I. RESULTS 


The HEA is a system of 2 k Ising spins Si = ±1, with 
Hamiltonian 

2 k -l 2 k -l 

H k [S] = - - Z h ^ (!) 

i<j =0 2=0 

where {J,; ; } are independent and identically distributed 
(IID) random variables defined as follows. Given two sites 
i,j, we denote by dij the hierarchical distance between i 
and j, i.e. the number of hierarchical levels that we need 
to ascend in the hierarchical tree starting from spins i, j to 
find a root common to i and j, see Fig. [TJr. In Eq. (HJ, we 
then chose Jij to be different from zero with probability 

Pij = 2~ 2 ^~ l \ (2) 

where nonzero Jij s are equal to ±1 with equal probabil¬ 
ity, the parameter a determines how fast spin-spin in¬ 
teractions decrease with distance, and we chose 1/2 < 
a < 1 because this is the interval where the hierarchical 
Edwards-Anderson model is expected to have a finite- 
temperature phase transition [19J- The magnetic fields 
{hi} in Eq. (JT]) are also IID random variables: hi is dif¬ 
ferent from zero with probability p k , and nonzero /i,s are 
equal to ±1 with equal probability. 


Since the probability @ with which non-zero bonds 
Jij are drawn in the Hamiltonian m is a function of the 
hierarchical distance dij, the correlation function [21J ] 

Ftf =m{SiSj) - (SMSj)) 2 ] (4) 

depends on i and j only through dij , where {) denotes 
the Boltzmann average with Hamiltonian ([!]) and inverse 
temperature /3 = 1/T. To obtain the correlation len gth , 
we recall that the hierarchical distance dij is related [22j 
to the dyadic norm | L by the relation 

\m-i{j)\ 2 = 2 k - d % ( 5 ) 

where 

k -1 

l{i) = J2ak-i-j2 j , i = 0,... ,2 k — 1, ( 6 ) 

i=o 

and the coefficients {aj} are given by the expression in 
base two of i: 

k -1 

* = E a j 2 7 ' ( 7 ) 

i=o 

Equation (O implies that the correlation function Tij de¬ 
pends on i, j only through | I(i) —I(j )| 2 : 

l{j)\ 2 ). (8) 


The choice @ for the probability of placing a bond im¬ 
plies that the average interaction strength on hierarchical 
level l > 1 between spins Si and Sj is 

E[4]=2- 2ct(i - 1 ), (3) 


We now relabel the lattice indexes i —» X(i) and we 
perform the Fourier transform of Tij @: this leads to the 
definition of the momentum-dependent spin-glass suscep¬ 
tibility 


where E[] denotes the average with respect to all ran¬ 
dom variables. Equation ([3]) shows that, for a given a, 
the power-law dependence on l of the average interac¬ 
tion strength is the same as in a HEA with Gaussian 
couplings studied recently [Q. The observation above 
suggests that the HEA introduced here and the HEA 
with Gaussian couplings have the same critical features: 
this will be the working hypothesis of the rest of the pa¬ 
per. In particular, this assumption implies that the two 
models share the existence, or the absence, of a spin- 
glass transition, and that they both have a MF behavior 


XsgO) 


^ ex P (^(* “ J')) T miU) ( 9 ) 

^ J2 ex p (^(* -■?')) /(i*— jI 2 )» 

i,j =0 ^ ' 


where p = 0,..., 2 k — 1 is an integer momentum variable, 
and in the second line we used Eq. ((SJ) and the identity 

a m) = i- 
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To extract from Eq. an}, we consider the mo¬ 
mentum values with the smallest dyadic norm: these are 
p = 0 and p = 2 fc_1 , and their dyadic norms read |0| 2 = 0 
and |2 fc ” 1 | 2 = 2 1 ~ k respectively [22}. The correlation 
length is then obtained by considering Eq. (flUl) for 
these values of p: 


t _ 0 k- i ( Xsg(O) _ 1 

& Vxs G (2 fe - 1 ) 


( 11 ) 


b 

+1 


-1 



We will now discuss a convenient way of computing the 
correlation length in terms of products of spin overlaps 
[24} , and thus extract from the MC data. We introduce 
a set of replicas {5'“} of the spin configuration S, and we 
set 


2 fc _ 1 

Qab(p) S i S i eX P T -• 

2=0 ' ' 

By using Eqs. ©, 02}, we have 


( 12 ) 


FIG. 1. Hierarchical Edwards-Anderson model in a magnetic 
field (HEA) with k = 3. (a) Structure of spin-spin interac¬ 
tions: each dot represents a spin, and arcs represent interac¬ 
tions between spins below them. The green path goes from 
the bottom to the top of the hierarchical tree until a common 
arc between the two circled spins is found: since this requires 
ascending two hierarchical levels, the hierarchical distance be¬ 
tween the two spins circled in green is dij = 2. (b) Long- 
wavelength Fourier modes: real part of the Fourier modes 
exp(27ripl(i)/2 fe ) for the HEA (black solid lines), where the 
index i = 0,... , 2 fc — 1 runs over the lattice sites from left to 
right, and the Fourier modes take values 1 and —1. We also 
plot the real part of the Fourier modes exp(27ripi/2 fc ) for a 
one-dimensional spin glass where spin-spin interactions decay 
according to the Euclidean distance | i — j\ (red dashed lines). 
Top: first Fourier mode, corresponding to p = 0 for both 
the HEA and for the power-law spin glass. Bottom: second 
Fourier mode, corresponding to p = 2 fe_1 for the HEA and to 
p = 1 for the power-law spin glass. 


XSG (p) = 2 fc {E {(QMp)QUp))} - 2E KQMp)QUp))} + 
+E KQMQUp))}}, (13) 


where in what follows the average () of a function of 
S 1 , S 2 , ... denotes the Boltzmann average over all repli¬ 
cas [241 1 - each replica having an independent Boltzmann 
measure with Hamiltonian © and inverse temperature 
/?. We now introduce the overlaps in the left and right 
half of the lattice 


2 —1 


2—1 


Qab^^Z T E Qab^^Z I E S i S - 


2 k ~ 


i=0 i—2 k ~ 1 

and by using Eqs. 03, 03, we have 


(14) 


1 2 - 1 nL , 

Qab( 0) =Q ab = ¥ J2 S ? S i = 

2=0 


Q R 


ab 


(15) 


In spin-glass models with nearest-neighbor interac¬ 
tions, the correlation length is extracted [23] from the 
Fourier transform of the correlation function by consid¬ 
ering the values of p with the smallest norm |p|. For the 
HEA, the natural definition of distance between sites on 
the lattice is not the Euclidean norm |p|, but the dyadic 
norm |p| 2 0: to extract the long-wavelength behavior, 
it is thus natural to consider the values of p with the 
smallest value of |p| 2 . To this end, we observe that the 
spin-glass susceptibility © is the Fourier transform of a 
function / which depends on i — j only through its dyadic 
norm, thus xsg(p) depends on p only through the dyadic 
norm |p| 9 of the momentum [22} . The long-wavelength 
behavior is then encoded into the limit |p| 2 —> 0, in which 
Eq. © takes the modified Ornstein-Zernike form [21] 


XSG (p) oc 


1 


>.-(2(7-1) 

s k 


I |2cr—1 

bl 2 


( 10 ) 


Qab( 2 fe 2 ) = 7 ^ E S?^exp(7ril(f)) 


(16) 


2=0 
2 fe — 1 


= ¥ E s?s!(-i) x «- 


2=0 

^2 fc_1 — 1 


2 K —1 


_ 1 { \ ^ qa qb \ A qa qb 

— 9 k I 2^ 2-^ 

\ 2=0 

Qa b ~Qab 


i=2 k ~ 1 


where in the third line we observed that I(i) is even if 
0 < i < 2 fc_1 - 1 and odd if 2 fc " 1 < i < 2 k - 1, see Eqs. 
©, ©. By using Eqs. 03, 03}, 05} , 03} we obtain 
the final expression for the correlation length 
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& = 2 fc_1 


/ E[{(Qf 2 + Qf 2 ) 2 )] - 2E[((Q& + Qf 2 )(Q& + Q&))] + n((Qr 2 + QfzXQk + Q&))] 
U[((Qf 2 ^ Qf 2 ) 2 )] - 2E[((Qf 2 - Qf 2 )(Q^ 3 - <?£)>] +E[((Qfa - <?£)(<& - Q&)>] 


(17) 


The derivation above provides a numerically convenient 
way of extracting the correlation length: to compute 
we do not need the full set of spin-spin correlations {T^}, 
but only the products of the left and right overlaps in Eq. 

ED- 

The first two long-wavelength Fourier modes which de¬ 
termine i.e. exp(27ripl(*)/2 fc ) for p = 0,2 fc 1 , are 
depicted in Fig. [TJx in particular, the second Fourier 
mode is (— 1) X W, and it is equal to 1 and —1 in the left 
and right half of the lattice respectively, see Eq. ED- 
This particular form of the Fourier modes is due to the 
hierarchical structure of interactions in the model: for 
example, a one-dimensional spin glass where interactions 
depend on the Euclidean distance rather than on the hi¬ 
erarchical distance (fj possesses a different second Fourier 
mode, given by exp(27ri i/2 k ). 


B. Numerical results for the correlation length and 
overlap probability distribution 


We will now study the existence of a phase transition 
by analyzing the behavior of the correlation length as a 
function of temperature: namely, we will assume that if 
there exists a critical point, then for T close to the critical 
temperature T c and large en oug h k the following finite- 
size-scaling relations hold [!). Toj 

a/2 fc "/ 3 = <?MF(2 fe/3 (T - T c )), 1/2 < a < 2/3, (18) 
£ k/2 k = 5 NM F (2 fc/l/ (T - T c )), 2/3 < a < 1, (19) 

where the correlation-length critical exponent exponent 
v is equal to 1/(2<r — 1) for 1/2 < a < 2/3 0 and ^mf, 
Snmf are two different scaling functions in the MF and 
NMF region respectively. 

We computed the correlation length ED with MC sim¬ 
ulations performed with the parallel-tempering algorithm 
l2f} . For zero magnetic field the Hamiltonian m is sym¬ 
metric under a global spin flip, thus the second and third 
term in both the numerator and denominator of Eq. ED 
vanish and can be computed by simulating only two 
replicas per temperature, while for nonzero field four 
replicas per temperature are needed. We combined par¬ 
allel tempering with the asynchronous multispin-coding 
method, where 64 disorder samples are simulated simul¬ 
taneously by encoding the values of the couplings Jij into 
a 64-bit integer [26}. This method allowed for equilibrat¬ 
ing an extensive number of samples 8 x 10 3 < S < 16 x 10 3 
and system sizes 2 6 ,..., 2 13 . 

Let us first focus on the zero-field case ph = 0: in the 
left panel of Fig. [5] we show £fc/2 fcw / 3 vs. T in the MF 
region, a = 0.6, while in the right panel we plot £fc/2 fc 


vs. T in the NMF region, a = 0.68. In the left and 
right insets we plot £&/ 2 kv / 3 and £*,/ 2 k respectively vs. 
the logarithm of the number of MC sweeps, in order to 
provide an equilibration test for the data in the main 
panels. The correlation-length curves cross at a finite 
critical temperature, showing that there is a spin-glass 
transition both in the MF and in the NMF case. This 
occurrence of a spin-glass transition in zero field is in 
agreement with a previous RG analysis for a HEA with 
Gaussian couplings and zero magnetic field, where a crit¬ 
ical FP has been shown to exist both in the whole MF 
region 1/2 < a < 2/3 and in the NMF region a = 2/3 -F e 
by means of an e-expansion [7]. 

To study the nonzero-field case, we need to choose a 
value for the probability ph of placing a local magnetic 
field: it can be shown that a suitable value of ph is ob¬ 
tained as a tradeoff between two effects. First, for large 
Ph the critical temperature—if any—is too low to achieve 
equilibration for sufficiently large system sizes. Second, 
for small ph a phase transition might appear, but this may 
be a spurious effect of the transition for ph = 0 shown in 
Fig. [2] rather than an actual transition in a finite field. 
In this regard, a natural way of checking whether the ap¬ 
parent transition is an artifact of zero-field effects is to 
compute the overlap distribution 

p(q)=E((S(q^Q 12 ))}. (20) 

If there is a phase transition, in the thermodynamic limit 
p(q) must converge to S(q — go) at the critical point [24} . 
where go = E[(Qi 2 )], and go is positive because 

® = iE E [^)] = ^EE[ ( m 

•£—0 2—0 

When the system size is finite, the effect of the zero-field 
transition arises in the shape of a long exponential tail 
of p(q) for q < 0 M: thus, a natural way of ensuring 
that we are in the true nonzero-field regime is to choose ph 
large enough that the support of p(q) is mostly localized 
for positive q, i.e. a q /q 0 < 1, where = E[((Qi 2 -g 0 ) 2 )]. 

We have thus chosen the value of ph according to the 
tradeoff above. In the MF region, a = 0.6, we take ph = 
0.36: the left panel in Fig. [3] shows that for this value of 
Ph the curves of £*,/ 2 kv ^ vs. T cross at a finite critical 
temperature T c . In addition, the right panel of Fig. [3] 
shows that for this value of ph the zero-field effects are 
weak: According to the top inset, the overlap distribution 
at T c has the typical left exponential tail due to zero- 
field effects [ 23 , nevertheless the main panel illustrates 
that the support of p(q) is mostly localized at positive q, 
and the bottom inset shows that <r q /qo <C 1 for large k. 
The same qualitative scenario arises in the NMF region, 
<t = 0.68, see Fig. [4] The left panel indicates that the 
curves £*,/ 2 k cross at a finite critical temperature, hinting 
at the existence of a phase transition in a field, and the 
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FIG. 2. Correlation length for zero magnetic field, ph = 0. Left panel: ^k/2 kv ^ 3 vs. temperature T in the mean-field 
region, a = 0.6 < 2/3, for system sizes 2 , with k = 6,..., 13 and v = 1/(2<t — 1): the curves cross at a critical temperature 
T c = 2.38 ± 0.24, which is also marked. Inset: £k/‘2 kl/ ^ 3 vs. log 2 t for the lowest temperature shown for each system size in the 
left panel, where t is the number of Monte Carlo sweeps, and for every t we plot £*,/2 fel '/ 3 computed from the last t/2 sweeps. 
Right panel: £k/2 k vs. T in the non-mean-field region, a = 0.68 > 2/3, with a critical-temperature value T c = 1.65 ± 0.11. 
Inset: £,k/2 k vs. log 2 t for the lowest temperature shown for each system size in the right panel, where £ k/2 k is computed as in 
the left-panel inset. 


T c 




FIG. 3. Correlation length and overlap probability distribution for nonzero magnetic field in the mean-field region, a = 0.6 < 
2/3, with ph = 0.36. Left panel: £*,/ 2 kl/ ^ 3 vs. temperature T and system sizes 2 fc , with k = 8,..., 13. Smaller sizes are affected 
by strong finite-size effects, and they are not shown in order to display clearly the curve crossings for larger sizes. The curves 
cross at a critical temperature T c = 1.46 ± 0.07, which is also marked. Inset: £&/ 2 kv ^ vs. log 2 t for the lowest temperature 
shown for each system size in the left panel, where t is the number of Monte Carlo sweeps, and for every t we plot £*,/ 
computed from the last t/2 sweeps. Right panel: p(q) vs. q for T = T c and same system sizes as in the left panel. Top inset: 
same plot as in the right panel in logarithmic scale on the vertical axis. Bottom inset: ratio between the standard deviation a q 
and the mean qo of p(q) as a function of k. 


right panel illustrates that such a transition is not an 
artifact due to the effects of the zero-field transition. 


C. Hallmark of non-perturbative effects 


The correlation-curve crossings in the right panel of 
Fig. U suggest that there is a spin-glass transition in 
the NMF region a > 2/3: interestingly, this result is 
at variance with a recent perturbative RG analysis for a 
HEA with Gaussian interactions [l8j, which shows with 


an e-expansion [I] that for a = 2/3 + e there is no per¬ 
turbative FP corresponding to a spin-glass transition in 
a field. The observation above thus raises the possibil¬ 
ity that the phase transition in the right panel of Fig. Q] 
results from non-perturbative effects, a scenario that we 
will investigate further in what follows. 

The RG study mentioned above shows that if there is 
a spin-glass transition described by a perturbative FP, 
then at the critical point the rescaled spin-glass suscep¬ 
tibility xsg 2 k( ' 1 ~ 2<T ' > must con verg e to a finite value for 
k —>• oo, where xsg = Xsg(0) [13]. The same RG anal¬ 
ysis demonstrates that in the MF region there is a crit- 
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FIG. 4. Correlation length and overlap probability distribution for nonzero magnetic field in the non-mean-field region, a = 
0.68 > 2/3, with ph = 0.16. Left panel: £i,/2 fc vs. temperature T and system sizes 2 k , k = 9,..., 13. Smaller sizes are affected 
by strong finite-size effects, and they are not shown in order to display clearly the curve crossings for larger sizes. The curves 
cross at a critical temperature T c = 0.88 ± 0.07, which is also marked. Inset: £i,/2 fc vs. log 2 t for the lowest temperature shown 
for each system size in the left panel, where t is the number of Monte Carlo sweeps, and for every t we plot £k/2 k computed 
from the last t/2 sweeps. Right panel: p(q) vs. q for T = T c and same system sizes as in the left panel. Top inset: same plot 
as in the right panel in logarithmic scale on the vertical axis. Bottom inset: ratio between the standard deviation a q and the 
mean qo of p(q) as a function of k. 


T? 



T? 



FIG. 5. Rescaled spin-glass susceptibility as a function of temperature for nonzero magnetic field. Left: xsg 2 k ( 1 ~ ~ 2 0 vs. T 
in the mean-field region, a = 0.6 < 2/3, for ph = 0.36 and system sizes 2 k , with k = 6,..., 13. Equilibration tests for xsg 
have been performed along the same lines as in Figs. [2] [3] [4] and they are not shown. The curves cross at a finite temperature 
T* = 1.71 ± 0.06, which is also marked. Right: same plot as in the left panel in the non-mean-held region a = 0.68 > 2/3 for 
p h = 0.16 and k = 6,..., 13. The crossing temperature of two curves with system sizes 2 k and 2 k+1 is an increasing function of 
k, and the estimated infinite-volume critical temperature is T* = 0.88 ± 0.06. 


ical FP, hence in this region the rescaled susceptibility 
XSG 2 fe ( 1_2cr ) must converge to a finite value at the critical 
temperature. This scenario is confirmed by the numerical 
data in the left panel of Fig. [5j which shows that in the 
MF region xsg 2 fe ( 1_2cr ) converges to a finite FP at a finite 
temperature T*-, which is within ~ 17% from the critical 
temperature T c obtained from the correlation length in 
Fig. [3] The situation is completely different in the NMF 
region: there, the RG approach shows that there exists no 
physical FP, thus suggesting that xsg 2 fc ( 1 ~ 2cr ) does not 
have a finite FP at any temperature [18|. Still, the right 
panel of Fig. [5] illustrates that XSG 2 fe ( 1-2<T ) converges to 
a finite FP at a finite temperature T*. In addition, T'*■ 


coincides within error bars with the critical temperature 
T c obtained from the correlation length in Fig. [I] 

Overall, the discrepancies between the predictions of 
the perturbative RG analysis and the numerical results 
in the NMF region raise the possibility that the phase 
transition resulting from the numerics in the NMF re¬ 
gion is associated with a non-perturbative FP. In this 
regard, a general feature characterizing a critical FP is 
the exponent v governing the exponential departure of 
RG trajectories as the temperature is moved away from 
its critical value [l]. To study the existence of a non- 
perturbative FP, we estimate v by means of the finite- 
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2*7 3 (T - T c ) 


2 k / v (T-T c ) 


FIG. 6. Estimate of the critical exponent v with a nonzero magnetic field. Left: £ y k/2 kl '^ 3 vs. 2 k / 3 (T — T c ) in the mean-field 
region, a = 0.6 < 2/3, for system sizes 2 k , k = 6,..., 13, where v = 5.0 and T c is determined by Fig. [3] The gray line envelopes 
the curves for largest system sizes, and it provides an estimate of the finite-size-scaling function <?mf. The value of v has been 
chosen to achieve the best possible overlap between the curves enclosed in the gray envelope, yielding the estimate v = 5.0 ±0.4. 
Right: £j./ 2 k vs. 2 k / v (T — T c ) in the non-mean-field region, a = 0.68 > 2/3, for system sizes 2 k , k — 6,..., 13, where v = 5.0 
and T c has been determined by Fig. [4] The the finite-size-scaling function jnmf is also shown, and v has been estimated in the 
same way as in the left panel, providing v = 5.0 ± 1.0. 


size-scaling equations m , ©: we chose v in such a 
way that for large k the curves ^fe/2 fcy / 3 vs. 2 k / 3 (T — T c ) 
and £fc/ 2 k vs. 2 k / v (T — T c ) in the MF and NMF region 
respectively merge into a single curve. In the MF region, 
a = 0.6, we have v = 5.0 ± 0.4, see Fig. [6] left, in agree¬ 
ment with the RG classical value v = 1/(2 a — 1) = 5 
associated with the MF FP [ij. The situation is again 
very different in the NMF region: if we use the MF FP 
prediction v = 1 /(2cr — 1) to estimate v for a = 0.68, we 
obtain v ~ 2.78. This values differs by nearly a factor 
of two from the MC estimate v = 5.0 ± 1.0 in the right 
panel of Fig. [6j even if a is within only ~ 2% from its 
upper-critical-dimension value a = 2/3. This discrepancy 
suggests that the FP giving rise to the phase transition 
in the numerics cannot be regarded as a slight perturba¬ 
tion of a MF FP with classical exponents. Importantly, 
this last observation may be at the bottom of the overall 
disagreement between the numerics and the RG analy¬ 
sis: the advance here is that the RG approach may fail 
to capture the critical behavior in the numerics because 
it is based on the hypothesis that the FP governing the 
transition in the NMF region is a perturbation of the MF 
one. 


II. CONCLUSIONS 


We investigated the existence of non-perturbative ef¬ 
fects in spin glasses with a Monte Carlo (MC) study for 
an Ising spin glass in a field where pairwise couplings are 
disposed in a hierarchical way—the hierarchical Edwards- 
Anderson model in a magnetic field (HEA). The spin- 
interaction decay in the HEA is controlled by the pa¬ 
rameter er, which is reminiscent of the dimension d for a 
ferromagnetic Ising model [28|: in particular, the mean- 
field (MF) and non-mean-field (NMF) regions a < 2/3 


and a > 2/2, correspond to d > 4 and d < 4 respectively 
iT9| . The HEA is suitable for assessing the existence of 
non-perturbative effects for two reasons: First, the binary 
nature of couplings and magnetic fields makes the HEA fit 
for non-perturbative numerical approaches, such as MC 
simulations 0] . Second, the hierarchical structure of the 
couplings reproduces the same average interaction struc¬ 
ture as in a HEA with Gaussian couplings recently stud¬ 
ied with perturbative renormalization-group (RG) meth¬ 
ods [l8[. This RG analysis provided the complete set of 
perturbative RG fixed points (FPs): in the MF region, it 
was shown that there is a FP corresponding to a phase 
transition, while no such FP exists in the NMF region. 
In the MF region the MC results are in agreement with 
the RG predictions: the numerics show the existence of a 
finite-temperature transition where the correlation length 
and the spin-glass susceptibility become singular at large 
sizes, providing an estimate of the exponent v in accor¬ 
dance with its classical value [7J. Conversely, the MC 
and the RG approach markedly disagree in the NMF re¬ 
gion, where the MC data for both the correlation length 
and spin-glass susceptibility support the existence of a 
finite-temperature transition. Also, the MC estimate of 
v in the NMF region differs by nearly a factor of two 
from its classical value, even if a is within only ~ 2% 
from the value a = 2/3 separating the MF and NMF 
regions. This picture is at variance with other systems 
where perturbation theory is well-behaved such as the d- 
dimensional Ising model, where the MC estimate of v in 
d = 3, i.e. within 25% from the upper critical dimension, 
differs by only ~ 25% from its classical value |2J. Taken 
together, these findings hint that the spin-glass transition 
resulting from the numerics in the NMF region cannot 
be captured with a perturbative framework, raising the 
possibility that this transition is related to a FP which 
cannot be described as a perturbation of the MF one. 

These findings raise a variety of possible scenarios con- 












cerning the structure of the low-temperature phase. On 
the one hand, the occurrence of a transition in a field 
is associated with the replica-symmetry-breaking (RSB) 
scenario, a MF picture characterized by the existence of 
exponentially many low-lying energy states (24j. Thus, 
the evidence provided here in favor of such a transition 
hints that the MF picture is accurate at least in some 
part of the NMF region. On the other hand, our findings 
show that the phase transition in the NMF region cannot 
be captured by an approach where corrections to the MF 
picture are considered as a perturbation. This fact raises 
the possibility that the physical features of this transition 
are markedly different from the MF ones: this possibility 
is reminiscent of scenarios alternative to the RSB one, 
such as the droplet picture [l3| or other scenarios inter¬ 
mediate between the RSB and the droplet one [29j. In 
this regard, the structure of the low-lying energy states 
could be directly analyzed by probing the low-energy ex¬ 
citations above the ground state fihl jdof and by charac¬ 
terizing, for example, their stiffness [3l( and surface di¬ 
mension [32{. A interesting possibility would then be to 
study how the resulting picture for the low-energy land¬ 


scape relates to the existence of non-perturbative effects. 
Taken together, the points above provide an interesting 
direction that is worth exploring in future studies. 
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